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In the recent experiment by Cren et al. [Phys. Rev. Lett. 102, 127005 (2009)], no hysteresis for 
vortex penetration and expulsion from the nano-island of Pb was observed. In the present paper, we 
argue that this effect can be associated with the thermoactivated surmounting of the surface barrier 
by a vortex. The typical entrance (exit) time is found analytically from the Fokker-Planck equation, 
written in the form suitable for the extreme vortex confinement. We show that this time is several 
orders of magnitude smaller than 1 second under the conditions of the experiment considered. Our 
results thus demonstrate a possibility for the thermal suppression of the surface barrier in nanosized 
low-T c superconductors. We also briefly discuss other recent experiments on vortices in related 
structures. 

PACS numbers: 74.25.Bt, 74.25. Uv, 74.78.Na 

I. INTRODUCTION 

Vortex nucleation in an infinite type-II superconductor becomes energetically favorable at lower critical field H c \. 
It is however clear that real superconductors are always bounded, so that vortices can appear in a sample, at low 
temperatures, only by penetration through surfaces. This process was for the first time considered by Bean and 
Livingston^ within the framework of the London approximation. In their approach, single vortex was represented by 
a straight line entering a semi-infinite sample through its surface. Bean and Livingston have found a surface barrier 
preventing such an entry. The barrier disappears when H is close to the thermodynamical critical field H c , which is 
much larger than H c \ in superconductors with high values of the Ginzburg-Landau parameter k (H c /H c i ~ k/Iiik,). 
Physically, Bean-Livingston barrier is a consequence of the competition between the two factors: (i) the interaction 
between the vortex and the Meissner current, which tries to push the vortex line inside the sample; (ii) the Magnus 
attraction of a vortex to the surface, which can be understood in terms of the interaction between the vortex and its 
mirror image ("antivortex"). Similarly, surface barrier exists for the vortex exit from the sample. Bean-Livingston 
barrier plays an important role in the vortex dynamics near the superconductor surface and leads to pronounced 
hysteresis effects^. 

As any potential barrier, Bean-Livingston barrier can be surmounted by thermal activation. This phenomenon 
was addressed by Petukhov and ChechetkirP, who have considered vortex nucleus as a half-loop spreading from 
the surface into the bulk. The profile of the surface barrier was calculated within the London approximation and 
the typical penetration time was estimated by using Fokker-Planck equation. The resulting formula is of Arrhenius 
type with the height of the surface barrier entering through the exponent. After substituting parameters typical for 
conventional low-T c materials, this height was estimated to be as large as 10 3 -10 5 in terms of fc^T, from which it was 
concluded that thermally activated vortex penetration is practically unobservable. 

However, two decades later, after the discovery of high-T c materials, thermal activation over the surface barrier has 
attracted great attention. It was showiPby Kopylov et al. that thermal surmounting of the barrier becomes possible 
in layered high-T c superconductors, where the mixed state structure is represented by quasi-2D pancake vortices. 
Burlachkov^ has analyzed the same phenomenon for high-T c superconductors of YBaCuO type with 3D vortex lines. 
He has demonstrated that the height of the surface barrier in these compounds is dramatically suppressed, compared 
to low-T c materials, to 10 — lOOfcsT, which leads to experimentally observable rates for thermally activated vortex 
penetration and expulsion. Indeed, these conclusions agree with numerous experiments, see e.g. RePj, as well as the 
theoretical papeiP 

Thus, it is generally believed that thermally activated flux penetration and expulsion in low-T c superconductors 
is not possible, while for the case of high-T c superconductors it is quite standard. However, in the very recent 
experiment!^ by Cren et al. no hysteresis for vortex entry and exit was found in a superconducting nano-island of Pb 
having an extremely small thickness (« 5.5 nm). Lateral dimension of this sample was so small (~ 100 nm) that it 
could accommodate only one vortex before turning to the normal state. This regime was referred in RefP as ultimate 
vortex confinement. 

The aim of the present paper is to show that thermal fluctuations can be responsible for the effective suppression 
of the surface barrier observed in the experiment of CrerP et al.. This cannot be firmly proved by using previously 
developed approaches to the thermal activation over the surface barrier, since all these schemes are based on the 
London theory, while this theory is not applicable for small-sized samples with lateral dimensions of the order of the 
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vortex core size. Moreover, under these conditions, vortex position is not an appropriate variable, because an initial 
stage of vortex nucleation is associated with a "nascent" vortex, considered for the first time by Kramer^. "Nascent" 
vortex represents an area at the sample's edge, where the order parameter is suppressed, but it does not vanish yet. 
As we show, the height of the potential barrier, under the experimental conditions of RefP, is mostly due to the 
vortex core nucleation rather than because of the motion of already formed vortex towards the island center. Thus, 
the origin of the surface barrier in the limit of a very small sample size is not identical to the one for the traditional 
Bean-Livingston barrier. This feature has already been discussed in RefM 

In the present paper, in order to circumvent these difficulties, we use the following alternative scheme, (i) We first 
solve Ginzburg-Landau equations approximately by using lowest Landau level approximation for the order parameter 
and determine the surface barrier in terms of the populations of corresponding levels. The island we consider is 
so small that it can accommodate only one vortex. Hence, to describe a barrier with a reasonable accuracy, it is 
sufficient to take into account only two levels (with winding numbers and 1), which correspond to the vortex-free 
and one- vortex states. Populations of these levels are then used as "good" variables instead of the vortex position. 
Notice that there are no essential difficulties in extending the number of variables by taking into account more levels, 
which can be necessary when considering larger structures, (ii) We then estimate "viscosity" coefficients associated 
with the motion of order parameter projected on the two Landau levels. Demagnetization effects are included into 
the consideration, (iii) At the final step, we solve the backward Fokker-Planck equation, again in terms of the Landau 
levels populations. Our main result is an analytical expression for the typical time for thermally activated vortex entry 
(exit). Notice that the lowest Landau level representation for the order parameter in mesoscopic superconductors was 
proposed in RefsP^^ and then used very widely. 

Our calculations show that the height of the surface barrier, under the conditions of the experiment considered, is 
approximately 20kgT that is more typical for high-T c materials than for conventional low-T c superconductors. The 
major reasons for the potential barrier to be not so high are (i) very small island thickness, (ii) very small island 
lateral dimensions which result in the substantial suppression of the order parameter, (iii) relatively high sample's 
temperature^ (~ 4.3 K). The preexponential "attempt time" is quite small - mainly because of the small system sizes, 
which result in a weak magnetic response of the island. This time was estimated to be of the order of 10~ 15 — 10~ 14 s 
which is much smaller than the similar quantity- 4 ^ for high-T c layered superconductors with pancake vortices (~ 10 -12 
s). The resulting first passage time r both for the vortex entry and exit is also quite small, r ^10~ 5 — 1CP 4 s. 

Thus, the results we here present demonstrate that the thermally-assisted suppression of the surface barrier can 
occur in nanosized samples made of low-T c superconductors, which are already available for experimental investigations. 

Since nanosized superconductors might be perspective for technological application (see e. g. Re f 14 * 15 1 ) and are 
even already used (for instance, for single-photon detection purposes^), the control of fluctuations in such structures 
is quite important. 

This paper is organized as follows. In Section II, we study surface barrier within the Ginzburg-Landau theory 
In Section III we estimate viscosity coefficients, with taking into account demagnetization effects. In Section IV, we 
solve the backward Fokker-Planck equation and determine first passage times both for the vortex penetration and 
expulsion. We conclude in Section V. 

II. SURFACE BARRIER 

In this section, we determine approximately the height of the potential barrier for vortex exit and entry under 
the conditions of experiment. Our results are applicable as long as the sample's radius remains comparable to the 
coherence length £(T). 

A. Main parameters 

The superconducting island studied in RefP had a " smoothed" hexagonal shape with a small thinner region in its 
center, as seen from Fig. 1(b) of RefP. The smallest lateral size of the sample was 110 nm, while the thickness was 
d = 5.5 nm. Sample temperature was T = 4.3 K. The transition between the vortex-free and one-vortex states, as 
obtained by STS measurements, occurred at the applied field of 0.235 T. The transition was reversible within 3 %. 
The reversibility implies that at this field energies of vortex-free and one- vortex phases become equal. 

There was a disordered layer between the island and the Si substrate, which limitates^ the quasiparticle mean free 
part to I w 2d. The zero-temperature coherence length was estimated^ by using an expression for dirty superconductors 
~ a/Co^j where £o = 80 nm is the coherence length in bulk Pb. This yields « 30 nm. The coherence length £(T = A.3K) 
should be larger, but it is not easy to find its precise value theoretically because of the small system sizes (and nontrivial 
shape), which makes the microscopic physics quite complex-^1. However, it was possible to extract the coherence length 
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directly from the measured zero-bias conductance profile in the vortex core, which yielded the value of 40 — 45 nnfl^. 
In this paper, we assume a value of 48 nm. For the penetration depth, we use the usual expression^ for dirty 

superconductors, A(T) ~ 0.615 ^o\ f^^jr i j where Ao — 40 nm is the penetration depth in bulk Pb, T c = 7.2 K. This 

gives A(4.3if) ~ 102 nm. Notice that the use of the prefactor 0.615 is, in fact, rather relative. The Ginzburg-Landau 
parameter, corresponding to the above numbers, is k ss 2. When considering the diamagnetic response of the island, 
one has to keep in mind that effective K is much higheiPSl since d -C A(T). 

Although the island has a rounded hexagonal shape, it is actually quite similar to a disc (see Fig. 1(b) of Rcf. 8 ). 
Therefore, we model it by the disc. As for its diameter, we use the value of 150 nm, which corresponds approximately 
to the largest ("bottom") lateral dimension of the island. This value together with 48 nm for the coherence length 
will allow us to reproduce the experimentally found field for the transition between the vortex-free and single- vortex 
states (0.235 T). 

Note that the radius of the disc in terms of the coherence length is approximately 1.55. This regime indeed 
corresponds to the "ultimate vortex confinement": the disc is able to accommodate only one vortex before the 
transition to the normal state, see e.g. RefPHl. 



B. Basic estimates 



Let us make some basic estimates. We compare the thermal energy kgT with the energy required to create a vortex 
nucleus, i.e., to suppress the order parameter to zero within the volume ~ £(T) 2 d. The similar idea is applied when 
calculating the Ginzburg-Levanyuk numbeJ^. Within the Ginzburg-Landau theory, the energy gain due to the vortex 
nucleus is 

- H f^arfd, (i) 

2^o 

where H C {T) is the thermodynamical critical field, given by 

' 2ttV2£{T)\(T)' 1 ' 

For the experiment performed at T = we obtain the following estimate for the ratio: 

keT ( H ^T m2d ) 1 ~ 1(r2 ' (3) 

which is basically not so small. Notice that the order parameter is strongly suppressed in the vicinity of the transition 
between vortex-free and one-vortex states. This should reduce condensation energy and enhance the ratio (3) in 
several times. Also, we expect that the barrier height to be much lower than the above estimate due to an additional 
reason: vortex nucleation is associated mainly with the redistribution of the density of superconducting electrons 
inside the island rather than with the change of their total number. Thus, the ratio of fc^T and the barrier height 
has to be of the same order as in high-T c superconductors^, where rates for thermally activated flux penetration and 
expulsion are observable. Below we show that this is indeed the case. 



C. Calculation of the surface barrier 



We start with the dimensionless Ginzburg-Landau functional for the energy of the disc in the superconducting 
phase compared to the normal state: 

F = ^^i{TfdJJ rdr (" I/I' + \ I/I" + IH V - a) jf) , (4) 

where integration is performed over the disc cross section (cylindrical coordinate system is used), r — R/£(T) is the 
dimensionless radius of the disc, / is the dimensionless order parameter, a — h e r/2 is the vector potential, which has 
only an azimuthal component, h e = H / H C 2(T) is the external magnetic field measured in terms of the upper critical 
field H c2 (T) = <i>o/27r^(r) 2 , and all the distances are measured in units of £(T). Due to small disc sizes, the external 
magnetic field penetrates the sample almost totally, so that the magnetic field inside the disc is essentially equal to 
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the external field. Since the disc thickness d is much smaller than £(T = 4.3K), the order parameter does not vary in 
z direction. Hence, the problem both for the order parameter and barrier height is effectively two-dimensional. 
The order parameter can be represented as a Furrier expansion 

f(r,^) = J2c nV > n (r)e-^. (5) 

n 

In the vortex- free state, only one coefficient among c„'s is nonzero, namely cq. In the one- vortex state, a nonzero 
coefficient is c\. The disc we consider is so small (in terms of £(4.3-fQ) that it can accommodate only one vortex. 
Therefore, with good accuracy, we can map the order parameter on the subspace with only two nonzero components: 

f(r, ip) ~ c ip (r) + c 1 (f 1 (r)e~ lip . (6) 

The superposition (6) actually describes vortex exit or entrance. If we start with the vortex- free state (ci =0), an 
increase of c\ first leads to the suppression of the order parameter at some spot on the disc edge. This suppression can 
be interpreted as a formation of a vortex nucleus or a "nascent" vortex^. At certain value of ci, the order parameter 
in the center of this nucleus gets suppressed to zero, and the phase of the order parameter now changes by 2tt when 
turning around this topological defect. Thus, the real vortex is formed. Further increasing of c\ leads to the vortex 
displacement towards the disc center, with its position r v given simply by the condition: |co|<po( r «) — |ci|yi(r„), 
which follows directly from Eq. (6). If we now start with the one- vortex state (co = 0), an increase of Co results in 
the vortex displacement towards the disc edge until it leaves the island. 

Further simplifications come by using, for fo(r) and ipi(r), the eigen functions for the kinetic energy operator 
corresponding to winding numbers L = and L = 1, respectively, and to zero radial quantum number (nodclcss wave 
functions). The latter corresponds to the lowest eigen value for the kinetic energy operator, see e.g. Ref.^. This choise 
for ifo(r) and y>i(r) can be justified by noting that at the transition to the normal state, the first Ginzburg-Landau 
equation can be linearized, i.e., reduced to the eigen value equation for the kinetic energy operator (with the eigen 
value 1). Therefore, in the vicinity of the transition, <p n (r) can be closely approximated by the the eigen functions for 
the same operator. This regime is indeed realized in the superconducting disc under conditions of "ultimate vortex 
confinement" , where the transition between the one- vortex and vortex-free states occurs at high fields, H sa i/ c2 (T), so 
that the density of superconducting electrons is strongly suppressed. In cylindrical coordinates, eigen value equation 
reads 



1 d ( dtfL \ / h e r 



r- 



r dr \ dr J \ 2 

This equation can be solved analytically: 
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<pl hr _ t = £ l 1 pl- (7) 



l l/2 ( h e r 2 \ ( 'h e - e L h e r 2 . 

where $ is Kummer function. Eq. (7) has to be supplemented by the boundary condition at the superconduc- 
tor/vacuum interface 



d<fL 
dr 



= 0. (9) 



Values of El for given h e , ro, and L can be found numerically from Eqs. (8) and (9). 

We now treat Co and c\ as variational parameters and calculate disc energy, as given by Eq. (4), as a function of 
these two parameters. Thus the energy of the system is obtained as a biquadratic function of variational parameters 
Co and C\ . The resulting expression can be minimized with respect to Cq and c\ . Without loss of generality, Cq and 
ci can be taken real and positive. Although this method is not exact, we expect that it allows one to obtain rather 
accurate estimate of the potential barrier height for such a small island with the radius only slightly exceeding the 
vortex core size^. 

The energy of the island, obtained from Eqs. (4)-(9), is given by 

F = H ^ T ^((T) 2 d (c 2 a + c\ ai + 40o + 413! + 2 7C gc?) , (10) 
Mo 
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where (L — 0, 1) 



1"0 



a L = (s L - 1) / ip 2 L rdr, (11) 
Jo 

iplrdr, (12) 



7 = / VoVirdr. (13) 
Jo 

In the range of parameters corresponding to the ultimate vortex confinement, a^, (3l, 7 ~ 0.1 — 0.5. Let us stress 
that these numbers are internally dependent on the external field. 

The energy of the island in the vortex-free state F(L = 0) can be found from Eq. (10) by minimization of F with 
respect to Co, while keeping c\ = 0: 



F{L = Q) = ^MH i{ Tf d (^.\. (14) 
Mo V 4 Po / 

Similarly, the energy of the single-vortex state is given by 

F{L = l) = - H ^aTfd(§r). (15) 

Now, in order to find a height of the potential barrier for vortex entry, we start from the vortex-free solution, c\ = 0. 
We then begin to increase c\ and we minimize F with respect to Co at fixed c\. As a result, the following "optimal" 
dependence of F on c\ is obtained 
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Cl < 



V 



2/3i' 



(16) 



and similarly for the vortex exit, F(l —> 0). The height of the potential barrier, which is determined by the saddle 
point in (cq, c\) space is given by the minimum of F(0 — > 1) with respect to c±: 



Vm = ^(TTh^L. (17) 
» 4( £ -ft) 

A systematic study shows that the vortex-free state is unstable with respect to the vortex penetration at H ~ 0.255 
T; while the single- vortex state becomes unstable with respect to the vortex expulsion at H ~ 0.195 T. Thus, we see 
that the maximum hysteresis (in absence of fluctuations) is ~ 30 %. Note that the width of the hysteresis region 
expands rapidly with increasing R/£(T). 

In Fig. 1 we plot the dependence of the barrier height (both for vortex entry and exit), given in units of fc^T, 
on the external field H . We see that heights of both barriers are « 23 fc^T in the vicinity of H — 0.235 T, where 
energies of the two phases become equal. The ratio of barrier height and fc^T is of the same order as for high-T c 
materials. The major reasons for a relatively low height of the potential barrier, as compared to the thermal energy, 
are (i) very small island thickness and (ii) very small disc radius which results in the substantial suppression of the 
order parameter in the relevant range of applied fields. 

In Fig. 2 we show the calculated surface barrier for the vortex entry, given by Eqs. (16), as a function of C\. It is 
quite remarkable that almost the whole interval of c\ from c\ = up to the saddle point corresponds to the vortex core 
nucleation rather than to the motion of a formed vortex inside the island. Fig. 2 thus demonstrates that, under the 
experimental conditions of RefP, a vortex position cannot be an appropriate variable to analyze thermal activation 
over the surface barrier, due to very small island's lateral dimensions. 
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FIG. 1: The dependence of the surface barrier height on the external field H. 
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FIG. 2: Potential barrier preventing an increase of ci, i.e., a growth of the population of the Landau level with winding number 
1 at H — 0.235 T. This growth describes the formation of the vortex nucleus at the island edge with the subsequent motion of a 
formed vortex towards the island center. Dotted line separates the state with the "nascent" vortex, when the order parameter 
is still nonzero everywhere, and the phase with the real vortex. 

III. VISCOSITY 

When vortices move, local magnetic field inside a superconductor changes in time. According to the Faraday law, 
this change produces a nonzero electric field which leads to an energy dissipation. The dissipation can be interpreted 
in terms of a viscous force acting on moving vortices. 

The microscopic theory of this phenomenon is rather complicated. However, reasonable esti mate s for the viscosity 
force can be obtained by using phenomenological approach within the Ginzburg-Landau theorjfl^SM] j n this paper, 
in order to estimate viscosity associated with the motion of the order parameter, we apply similar ideas. 

When finding the island energy, we did not take into account that the magnetic field varies inside the sample. This 
is fully justified since the disc is very thin (d <C A(T)) and its radius is also small, R ~ A(T), that allows one to use a 
perturbation theory^ in terms of a small parameter l/re^y, where K e // is much larger than k due to demagnetization 
effects. Physically, this implies that the additional field generated by the Meissner current and by the vortex is small 
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in comparison with the external uniform field. If we know the solution of the first Ginzburg-Landau equation / for 
the order parameter in leading order, the vector potential a. add of this additional field can be found as 



_J_ f 
iadd ~4ir J 



-dV, (18) 
r 



where an integration is performed over the island volume and j is a supercurrent density given by 



J = 2 

K 



a|./f + g (/* V/ - /*V/) 



(19) 



where a = e v h e r/2 is the dimensionless vector potential in leading order which corresponds to the applied uniform 
field h e . Instead of using the exact solution for /, we take the lowest Landau level approximation, as done in Section 
II. 

Vortex entry or exit is associated with / changing in time. This change generates the local electric field e, which 
can be found from the Faraday law 

<9(rot & add ) 

rote= — , (20) 

and hence 

e = -^, (21) 

where a add is given by Eq. (18). 

The order parameter inside the island is already strongly suppressed in the vicinity of a transition to the one- vortex 
state. Therefore, the conductivity is close to the conductivity of Pb in the normal state <r n . Then, the dissipation 
rate in the sample W can be estimated as 

W/ (H c2 (T) 2 Z(T) 5 ) ~a n J e 2 dV, (22) 

where the denominator in the left-hand side comes from the normalization condition. 

Let us now consider the process of vortex entry within our representation for the order parameter, as given by Eq. 
(6). Vortex penetration is associated with the growth in time of c\ starting from 0, along the optimal path. It is then 
easy to see, from Eq. (19), that the leading-order in c\ contribution to j has the following components 

j<p = -^coci/i/o (h e r - ^ cos</?, 

Jr = — ^c ci (f[ r f - /o r /i)sin<p. (23) 

K 

We substitute Eqs. (23) to Eq. (18). Since we are interested in a add inside a thin disc, d <C R, we can neglect the 
dependence on z in the denominator of the integrand in the right-hand side (RHS) of Eq. (18). Then, for the radial 

component a^ d , after performing integration on angle, we obtain 

( r ) dcOS W , . ,„ 

a -=2<(^ C0Cl3(r) ' (24) 

where 



1 f r ° 

9(r)=~ rT 1 (/ lri /o-/^/i) 
' Jo 

where K and E are complete elliptic integrals of the first and second kind, respectively. Similarly, & a v dd is found, 
which is proportional to sin</?, while the integrand in the RHS of Eq. (25) now contains /o/i (h e r\ — l/^i) instead of 
(/in/o ~ /on/i)- The motivation to calculate viscosity only in leading order in c\ will be clarified in Section IV. 
By using Eq. (21), we find the local electric field e as a function of the time derivative c' lt 

er = -^ COS ^2^^\l-W 9{r) ' (26) 



-{n+r)E 



r + r\ 



+ 



2y/rn 



dn, 



(25) 



and similarly for e v . The dissipation rate is then determined from Eq. (22) as 

W~(c lt ) a n d 2 , (27) 

where 8 is found from integrating numerically the r-dependence of e, as given by Eq. (25), and by using calculated 
values of «o and Po- We then arrive to the simple estimate 

S~10- 2 . (28) 

The "force" of resistance against the growth of c\ is given by W divided by c' lt . This "force" is proportional to c' lt . 
The viscosity coefficient ri^\ i.e., the proportionality factor, reads 

(29) 

The same expression is found for the viscosity coefficient associated with the growth of Co (the case of a vortex exit). 
The estimate for S remains valid as well. 

We stress that rf Cl ' , by definition, is directly linked to the change of the population of the level with winding 
number 1 and not to the change of the vortex position. 

From Eq. (29) we see that proportional to d 3 and not to d, as can be naively expected. Such a strong 

dependence on d is due to demagnetization effects that significantly diminish the diamagnetic response of the sample 
and, hence, the electric field. This response is further suppressed due to extremely small lateral dimensions of the 
island, because the supercurrent created by the vortex is localized within a confined area ~ £(T) 2 . 

IV. FIRST PASSAGE TIME 

In Sections II and III we have found a profile of the potential barrier and the viscosity in terms of c\ (cq). We 
now are in a position to find the typical time for the vortex penetration (expulsion). This can be done by using the 
backward Fokker-Planck equation, also written in terms of C\ (c ). 

Let us first consider vortex entrance. The dependence U(a) of the system energy on c% is plotted in Fig. 2. 
Initially, the system is in the left-hand potential well (vortex-free state, a =0). To calculate the average time for the 
transition to the right-hand potential well (one-vortex state) we use the formalism developed in RefPP. 

We introduce the probability p(c[,t | Ci,0) to find the system at c[ in the moment of time t, provided it was at a 
in the moment t — 0. The relevant interval of possible values for a is between the two potential wells, as shown in 



Fig. 2 (ci <E [a, b], where a = 0, b = <J — ^-)< We also denote the saddle-point value of c\ as s. Then, the probability 
G(a , t) that the system is still in this interval is 

G(c x ,t)= [ Kc'i,* I ci,0)dc' v (30) 

J a 

Since the system is homogeneous in time, we can replace p{c[, t | Ci,0) by p(c' 1; | Ci, —t), for which the backward 
Fokker-Planck equation applies 

dp(c' 1 ,0\c 1 ,-t) = 1 dU dpje', , | ci, — f) k B T d 2 P (c'i 1 I ci, —t) 

d(-t) 7]^) da da 2^1) dc\ 

We can now switch back to p(c' 1: t | ci,0) and to integrate both sides of Eq. (31) on ci 6 [a,b]. Taking into account 
Eq. (30), we obtain an equation for G(a,t): 

dG(a,t) = l dU dG(g,t) k B T d 2 G(g,t) 

dt t/ci) da da 2ry( c i) def ' [ ' 

We also impose a condition 

p(c / l! 0|ci,0)=5(ci-ci) ) (33) 
which shows that at t = the system was in the definite state. From this condition and Eq. (30), we get 



G( C1 ,0)=1, 



(34) 
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which evidences that the system was inside the interval a <G [a, b] at t = 0. 

We then assume a reflecting condition at a = a, suggesting that the system never leaves [a, b] interval through this 
end. For the opposite limit, a = b, we impose an absorbing condition, which means that once the system reaches 
this point, it never enters again the [a, b] interval. We now calculate the average time r(ci), which corresponds to the 
time, system spends in the interval, provided it was at some fixed a in the initial moment. Thus, r(a) will give us 
the first passage time from the well at a = a (vortex- free state) to the well at a — b (one- vortex state). The time 
r(a) can be expressed from G(a,t). We first note that G(a,t) gives the probability that the system at the moment 
t is still within the interval [a, b], while G(ci,t + dt) is the probability that it is still there at t + dt. The difference, 
— dG< Qt'^ dt, represents the probability for the exit within the time interval [t, t + dt]. Then, the average exit time 
t(ci) is given by 

r(ci)=- f dt{^^)=- 1^ ^ dG( Cl ,t). (35) 



Using Eq. (32) for G(c\,t), we can obtain the differential equation for r(a) by integrating both sides of Eq. (32) 
on t from to oo: 

1 dU dr( Cl ) k B T 8 2 t(ci) _ , 



?y( c i) da da 2r/ c i) dej 

where we have taken into account Eq. (34) . Eq. (36) defines a first passage time from arbitrary a G [ a , b] to b. This 
equation must be supplemented by the conditions for the reflecting boundary at a and absorbing boundary at b 



dr(a) 



da 



Cl- 

T 



= 0, (37) 
(b) = 0. (38) 



The reflecting condition can be obtained from the vanishing of the probability current, while absorbing condition 
means that the system immediately exits the interval, once it is initially put at the absorbing edge. 

Ordinary differential equation (36) with boundary conditions (37) and (38) can be solved by standard methods. 
Finally, for r(a), which will be below denoted as r, we obtain 



— ! exp (S?)/> <c,)exp (^ 1 '" 1 -'" 



When integrating over y, the dominant contribution to the RHS of Eq. (39) is given by those values of y for which 
exp (U(y)/kBT) is largest, i.e., by the neighborhood of the saddle point s (see Fig. 2). At the same time, the internal 
integral is essentially constant in this area. It is therefore can be approximated by replacing y in the upper limit of 
integration by b. The RHS of Eq. (39) is then rewritten through a product of two independent integrals 



(M^H(i>''-(^)4 



The very high accuracy of such a transformation can be easily verified numerically. It is directly linked to the fact 
that the barrier height is large compared to k B T. 

The dominant contribution to the first integral in the RHS of Eq. (40) is provided by the neighborhood of s, while 
for the second integral it is given by the neighborhood of a. The latter is basically the reason why, in Section III, we 
have calculated the viscosity coefficient only in leading order in a- taking into account higher-order terms does 
not lead to any significant changes. We now expand U in the vicinity of s and a as 

U(x *s)kU(s) -B°A( x -8)\ (41) 
U(x sa a)^U(a) + ^ U "^ ( x - a ) 2 . (42) 
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It is easy to see that for the simple potential, given by Eq. (16), |f/"(s)| = 2\U"(a)\. After substituting Eqs. (41) 
and (42) to Eq. (40) and integrating, we obtain 

r=r exp(^%), (43) 



where U en — U(s) — U(a) is the barrier height, as given by Eq. (17). The "attempt time" tq is 



U"(a) 

By using Eq. (16), for U"(a) we get 



7o~ TTjrh ■ ( 44 ) 



U"{a) = 2 J^i(Tfd U - , (45) 



Mo V A) 

with [oti — ao7//?o) — 0.1. 

Thus, Eqs. (44) and (45) together with Eq. (29) for rf Cl ' yield the following estimate for tq 

T ~cr n /z ^. (46) 

We may assume that the conductivity is linear in T, so from the known o~ n of Pb at room temperature we can estimate 
a n at the conditions of the experimeniP as a n > 10 8 fl~ 1 m~ 1 . Then, tq is — 10~ 15 — 10~ 14 s. This value together with 
our result for the surface barrier height, U en « 23fcsT, gives r — 10~ 5 — 10~ 4 s. The very small value for To is due to 
the fact that the supercurrent is circulating within the tiny volume ~ £(T) 2 d, so that the magnetic response of the 
island is quite weak. Moreover, the order parameter is significantly suppressed. In addition, the distance from the 
island's edge to its center, which vortex has to pass, is very short. 

The above derivation was presented for the case of a vortex entry. The evaluation of the typical time for the vortex 
exit is quite similar. We have found that Eq. (46) remains applicable, which leads to similar order-of-magnitude 
estimates for r. 

Experimentally, it has also been observed that the magnetic field corresponding to the transition between one- vortex 
and vortex-free states was reversible within few percents^. This can be understood by noting that, when changing 
the magnetic field away from 0.235 T, barriers for exit and entrance become asymmetrical. This means that the 
equilibrium probability distribution becomes extremely strongly peaked within the deeper well, as follows directly 
from Fig. 2. Indeed, it is easy to see that the difference in heights of the barriers gives the difference in energies 
between the bottoms of the two wells, because (U(s) — U(a)) — (U(s) — U{b)) = U(b) — U{a). Since the absolute value 
of the difference in heights far exceeds fc^T already in the vicinity of 0.235 T (as seen from Fig. 2), the probability 
to find a system in a more shallow well is ~ exp(— \U(a) — U(b)\ /fcgT) << 1. Thus, the system gets stabilized in the 
deeper well. 

Note that To, as given by Eq. (46), is not explicitly dependent on T, in contrast to the similar quantity calculated 
in RefP within the London approximation. Technically, this is due to the fact that the barrier profile, as given by Eq. 
(16), is proportional to the second power of the relevant variable (with which r\ is associated), i.e., to c\ at c\ — > 0, 
while the barrier profile in RefP is proportional to the first power of vortex displacement when this displacement is 
small. Through the expansion (42) , this affects the second integral in the RHS of Eq. (40) and therefore leads to the 
mentioned difference. Physically, it might be attributed to the fact that vortex penetration processes for the cases 
of a half-infinite superconductor and the island in the ultimate vortex confinement regime are qualitatively different. 
Indeed, in the first case this process is associated with the motion of a "rigid" vortex line inside the superconductor, 
while in the second case it is represented by the formation of a "soft" vortex nucleus at the edge of the disc. In the 
first case, the barrier is due to the pronounced magnetic interaction with the superconductor's surface. In the second 
case, it is due to the condensation energy increase and the residual magnetic interactions. 

We also wish to notice that the effect of thermal fluctuations of the order parameter on vortex penetration and 
expulsion in mesoscopic superconductors was previously studied^ by Hernandes et al. by the numerical integration 
of the time-dependent Ginzburg-Landau equations. In these computations no thermal activation of vortices over 
the surface barrier was observed for parameters corresponding to samples made of low-T c superconductors, while 
this activation was shown to be possible for those samples, which are made of high-T c materials. This is consistent 
with previous experiments on Al discs^ and micron-size structures made of high-T c materiaPR We think that the 
difference between the results of the present paper and of RefP^ concerning low-T c samples can be attributed to the 
fact that in calculations of RefP^ the main parameters were chosen in such a way as to model experiments^, where 
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(i) discs of much larger thickness were studied compared to the island of RefP (one order of magnitude and more), 

(ii) the order parameter was not significantly suppressed due to larger discs radii, (iii) sample temperature was lower 
in one order of magnitude or so compared to that of RefP (T c of Al is at least several times lower than T c of Pb) . 

In the very recent experiment!^, strong current-induced fluctuation phenomena were observed in ultrathin (w 6 
nm) superconducting nanowires fabricated in a form of a meander and made of NbN. Lateral dimensions of sample's 
typical segment, in terms of £(T), were one order of magnitude larger than that for the island in the ultimate vortex 
confinement regime. Structures of this kind are used in single-photon detectors, while fluctuations are major source 
of dark-count events in such devices. Fluctuation phenomena were explained by a thermally-activated vortex entry, 
as well as by the unbinding of vortex-antivortex pairs. The height of the surface barrier, as calculated in RefP^, was 
at least one order of magnitude larger (in terms of fc^T) than the one for the superconducting island studied in the 
present paper. Nevertheless, vortex penetration events were detectable. 

Ref.GSO deals with the experimental investigation of superconducting islands of Pb, which were quite similar to the 
islands studied in RefP and in the present paper. The thickness of those islands was even smaller than that of 
RefP. However, measurements were performed at significantly lower temperature of 2.0 K. Probably, this is why 
some hysteresis for vortex penetration and expulsion was observed. Indeed, the activation time depends exponentially 
on the height of the potential barrier this height being much larger than fc^T. Nevertheless, it was mentioned in 
RefP^l that the width of the hysteresis region was noticably smaller than the theoretical values obtained from the 
Ginzburg-Landau theory. The results of the present paper evidence that thermal fluctuations might be responsible 
for an observed shrinkage of the hysteresis region. 

Let us finally mention a very recent experimental on one-atomic layer Pb films, where magnetic vortices were 
directly visualized. Thermal fluctuations must be strongest in such ultimately thin films. 

Since thermal fluctuations are more pronounced in low-dimensional systems, an understanding of vortex matter 
fluctuations in thin and small superconducting structures is of certain importance. It is howevere clear that more 
intensive experimental, as well as theoretical efforts are needed to reveal systematically those confinement and tem- 
perature regimes which favor thermal fluctuations in such nanostructures or, on the contrary, which enable one to 
safely avoid them. 

V. CONCLUSIONS 

Although it is generally believed that the surface barrier in low-T c superconductors is too high to be suppressed by 
thermal fluctuations, in this paper we have shown that such a scenario is possible in small and ultrathin supercon- 
ducting samples studied in modern experiments. In particular, we have demonstrated that an absence of hysteresis for 
vortex entry and exit in the nanosized island of Pb observed in the very recent experiment of RefP (and the shrinkage 
of hysteresis region reported in the similar experiment of Ref.^) can be explained in terms of a thermal activation of 
a vortex over the surface barrier. 

Lateral dimensions of these islands were so small that they could accommodate only one vortex before gradual 
transition to the normal state. For this reason, London approximation, which was previously applied to describe 
thermal surmounting of a surface barrier by vortices, is not applicable. 

We therefore have developed an alternative theoretical approach based on the lowest Landau level approximation 
for the Ginzburg-Landau order parameter, which is suitable for such extremely small structures. The surface barrier 
profile was calculated in terms of populations of the two relevant Landau levels being used throughout the paper as 
"good" variables instead of the "bad" one (vortex position). The viscosity coefficient, associated with the motion of 
the projections of the order parameter, was estimated with account of demagnetization effects. Finally, the expression 
for the typical time of thermally-assisted vortex entry (exit) was obtained from the Fokker-Planck equation, also 
written in terms of the populations of Landau levels. This expression is of Arrhenius type. 

We have found that the barrier height, under the conditions of experiment, is nearly 20 in terms of a thermal 
energy fc^T, which is close to the similar quantity for high-T c materials, where thermal suppression of a surface barrier 
is observable. Although the exponent of this ratio is still very large, the preexponential "attempt time" is quite small, 
so that the typical time both for the vortex exit and entrance is several orders of magnitude smaller than 1 s. 
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